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The cross-correlation matrix between equities comprises multiple interactions between 
traders with varying strategies and time horizons. In this paper, we use the Maximum 
i Overlap Discrete Wavelet Transform (MODWT) to calculate correlation matrices over 

different time scales and then explore the eigenvalue spectrum over sliding time windows. 
The dynamics of the eigenvalue spectrum at different times and scales provides insight 
into the interactions between the numerous constituents involved. 

Eigenvalue dynamics are examined for both medium and high-frequency equity re- 
' turns, with the associated correlation structure shown to be dependent on both time 

' and scale. Additionally, the Epps effect is established using this multivariate method 

, and analysed at longer scales than previously studied. A partition of the eigenvalue 

time-series demonstrates, at very short scales, the emergence of negative returns when 
C "3 ' the largest eigenvalue is greatest. Finally, a portfolio optimisation shows the importance 

of time-scale information in the context of risk management. 

o . 

1. Introduction 

In recent years, the equal-time cross-correlation matrix has been studied extensively 
for a variety of multivariate data sets across different disciplines such as financial 
data [H [H El SI El El El [H (HI [13 HI] , electroencephalographs (EEG) recordings QH 
Unim], magnetoencephalographic (MEG) recordings [T5] and others. In particular, 
Random Matrix Theory (RMT) has been applied to filter the relevant information 
from the statistical fluctuations inherent in empirical cross-correlation matrices, 
constructed for various types of financial, [2[1|31I1[51I1[71[E1G5], and MEG ' 
[16j . time-series. In this paper, we extend previous Complex Systems research, by 
examining time and scale dependent dynamics of the correlation matrix, in order 
to better understand the ever-changing level of synchronisation between financial 
time-series. 

The dynamics of the largest eigenvalue of a cross-correlation matrix over small 
windows of time was studied, 17J, for the Dow Jones Industrial Average (DJIA) and 
DAX indices. This analysis revealed evidence of time-dependence between 'draw- 
downs'/'drawups' and an increase/decrease in the largest eigenvalue. The inter- 
actions between stocks of two different markets (DAX and DJIA) were then in- 
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vestigated, |18j , revealing two distinct eigenvalues of the combined cross-correlation 
matrix, corresponding to each market. Adjusting for time-zone delays, the two eigen- 
values coincided, implying that one market leads the dynamics in the other. 

It has been suggested recently by several authors, [El [20l EU [22] [23] that there 
may, in fact, be some real correlation information hidden in the RMT defined 
part of the eigenvalue spectrum. A multivariate technique involving the equal-time 
cross-correlation matrix has been shown, [12] . to characterise dynamical changes 
in nonstationary multivariate time-series. Using this technique, the authors also 
demonstrated how correlation information can be observed in the noise band. As 
the synchronisation of k time-series within an M— dimensional multivariate time- 
series increases, a repulsion between eigenstates of the correlation matrix results, in 
which k levels participate. Through the use of artificially-created time-series with 
pre-defined correlation dynamics, it was demonstrated that there exist situations, 
where the relative change in eigenvalues from the lower edge of the spectrum is 
greater than that for the large eigenvalues, implying that information drawn from 
the smaller eigenvalues is highly relevant. 

This technique was subsequently applied to the dynamic analysis of the eigen- 
value spectrum of the equal-time cross-correlation matrix of multivariate Epileptic 
Seizure time-series, using sliding windows. The authors demonstrated that informa- 
tion about the correlation dynamics is visible in both the lower and upper eigen- 
states. A further study, [TS], which investigated temporal dynamics of focal onset 
epileptic seizure^, showed that the zero-lag correlations between multichannel EEG 
signals tend to decrease during the first half of a seizure and increase gradually be- 
fore the seizure ends. Information about cross correlations was also found in the 
RMT bulk of eigenvalues, [T3], with that extracted at the lower edge statistically 
more significant than that from the larger eigenvalues. Application of this technique 
to multichannel EEG data showed small eigenvalues to be more sensitive to detec- 
tion of subtle changes in the brain dynamics than the largest. A time and frequency 
based approach, has recently demonstrated the importance of correlation dynamics 
at high-frequencies during seizure activity, [14] . 

Although originally applied to complex EEG seizure time-series, [12] [13l [To] , 
these techniques were adapted to financial time-series, [11] . where the correlation 
dynamics for medium-frequency data were calculated using short time-windows. 
Changes in the correlation structure were shown to be visible in both large and 
small eigenvalue dynamics. Further, 'drawdowns' and 'drawups' were shown to be 
market dependent, with corresponding characteristic changes in relative eigenvalue 
size found at both ends of the eigenvalue spectrum. This analysis was carried out 
using daily data, with the effects of varying data granularity not examined. In this 
paper, we build upon this analysis by studying the multiscale effects on correlation 
dynamics using medium and high-frequency data. 



a A focal onset or partial seizure occurs when the discharge starts in one area of the brain and then 
spreads over other areas. 
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The correlation structure for various markets was studied in detail, |24j . for 
data encompassing a number of time horizons ranging from 5 to 255 minutes. By 
removing the centre of mass, (the market), the authors found that the correlation 
structure is well defined at high frequency, while for the original data the structure 
emerges at longer time horizons, (using stocks from the New York Stock Exchange). 

Wavelet multiscale analysis has been used by numerous authors to decompose 
economic and financial time-series into orthogonal time-scale components of varying 
granularities. Gengay and co-workers have variously examined the scaling properties 
of foreign exchange volatility, [25], and volatility models without intraday season- 
alities, |26j . using wavelet multiscaling techniques. The systematic risk in a Capital 
Asset Pricing Model was estimated over different granularities, [27], where it was 
shown that the return of a portfolio and its Betg(j become stronger as the scale 
increases for the S&P 500. 

An extensively studied characteristic of stock market behaviour, is the increase of 
stock return cross-correlations as the sampling time scale increases, a phenomenon 
known as the Epps effect, [25]. More recently, analysis of time-dependent correla- 
tions between high-frequency stocks demonstrated, however, that market reaction 
times have increased due to greater efficiency, [23] • A diminution of the Epps effect 
with time is one consequence of increased market efficiency. Trading asynchronicity 
was demonstrated to be not solely responsible for the effect, [3U], with the char- 
acteristic time apparently independent of the trading frequency. Further analysis 
using a toy model of Brownian motion and memory less renewal process, [31 , found 
an exact expression for the Epps frequency dependence, with reasonable fitting also 
for empirical data. In fact, the effect was shown, |32| . not to scale with market 
activity but to be due to reaction times, rather than market activity. 

In this paper, we extend the multivariate correlation technique, first applied 
to complex EEG seizure data, [T2J EH US] , and later to medium frequency finan- 
cial data, |llj . This interdisciplinary approach, allows characterisation of correlation 
changes in time. By decomposing equity returns into their component scales in short 
time- windows, we are able to study the correlation and associated eigenspectrum 
at various scales, allowing a time-scale analysis of the eigenspectrum. Using both 
medium and high-frequency stock returns, we examine the dynamics of the large 
eigenvalues and demonstrate the time-scale dependence of the correlation structure. 
We show how the Epps effect, [28], can be demonstrated through a multivariate ap- 
proach and expand previous work by examining considerably longer scales. Finally, 
we look at some applications of the work, with scale dependent risk characterisa- 
tion and portfolio optimisation examined. Sections ([2] — |3j) describe the techniques 
and data used, Section (f4| details the results obtained and conclusions are given in 
Section ©. 



b The Beta of an asset is a measure of the volatility, or systematic risk, of a security or a portfolio 
in comparison to the market as a whole. 
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2. Methods 

2.1. Wavelet multiscale analysis 

Wavelets provide an efficient means of studying the multiresolution properties of a 
signal, allowing decomposition into different time horizons or frequency components 
(Discrete Wavelet Transform, DWT), [331 131] • The definitions of the two basic 
wavelet functions, the father <f> and mother ip wavelets are given as: 



where j = 1, ... J in a J-level decomposition. The father wavelet integrates to 1 
and reconstructs the longest time-scale component of the series, while the mother 
wavelet integrates to and is used to describe the deviations from the trend. The 
wavelet representation of a discrete signal f(t) in L 2 (R) is: 



where J is the number of multiresolution levels (or scales) and k ranges from 1 to 
the number of coefficients in the specified level. The coefficients sj.k and dj tk are 
the smooth and detail component coefficients respectively and given by 



Each of the coefficient sets Sj, dj, dj-i, . . . d\ is called a crystal. 

In this paper, we apply the Maximum Overlap Discrete Wavelet Transform 
(MODWT), [33, 34 s , a linear filter that transforms a series into coefficients related 
to variations over a set of scales. Like the DWT it produces a set of time-dependent 
wavelet and scaling coefficients with basis vectors associated with a location t and a 
unitless scale tj = 2 J_1 for each decomposition level j — 1, . . . , Jq . The MODWT, 
unlike the DWT, has a high level of redundancy, however, is nonorthogonal and 
can handle any sample size N. It retains downsampledf) values at each level of the 
decomposition that would be discarded by the DWT. This reduces the tendency for 
larger errors at lower frequencies when calculating frequency dependent variance 
and correlations, fSection l2.3p . as more data is available. Having decomposed the 
equity returns into their component time-scales, we then calculate the correlation 
between these components. 



tf> jtk (t) = 2i <j> (2H - k) 
i/t jtk (t) =2*V(2*t-A) 



(1) 

(2) 




(3) 




(4) 




c Downsampling or decimation of the wavelet coefficients retains half of the number of coefficients 
that were retained at the previous scale and is applied in the DWT 
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2.2. Correlation dynamics 

The equal-time correlation matrix between time-series of stock returns is calculated 
using a sliding time window where the number of stocks, N, is smaller than the 
window size T. Given time-series of stock returns Ri (t), i = 1, . . . , N, we normalise 
the time-series within each window as follows: 



(t) = Mt^i (!) (6) 



where er.j is the standard deviation of stock i = 1, . . . , N and Ri is the time average 
of Ri over a time window of size T. Then, the equal time cross-correlation matrix, 
expressed in terms of r, (t) , is 

CtfSfaWr^t)) (7) 

The elements of are limited to the domain — 1 < Cy < 1, where = 1 defines 
perfect positive correlation, Cij = — 1 corresponds to perfect negative correlation 
and Cij — corresponds to no correlation. In matrix notation, the correlation is 
expressed as C = ^RR*, where R is an TV x T matrix with elements r.; t . 

The eigenvalues Xi and eigenvectors Vj of the correlation matrix C are found 
using Cvi = AjVj and then ordered according to size, such that Ai < A 2 < . . . < 
A at. Given that the sum of the diagonal elements of a matrix (the Trace) remains 
constant under a linear transformation, J^. A^ must always equal the trace of the 
original correlation matrix. Hence, if some eigenvalues increase then others must 
decrease, to compensate, and vice versa (Eigenvalue Repulsion). 

There are two limiting cases for the distribution of the eigenvalues [TH ITS] , 
with perfect correlation, Q w 1, when the largest is maximised with value N (all 
others taking value zero). When each time-series consists of random numbers with 
average correlation d ~ 0, the corresponding eigenvalues are distributed around 
1, (where deviations are due to spurious random correlations). Between these two 
extremes, the eigenvalues at the lower end of the spectrum can be much smaller 
than X max . By calculating the eigenspectrum associated with the correlation matrix 
at different scales in each time window, we get time-series of eigenvalues A mrl (r), 
where m denotes the eigenvalue number and n is the scale studied in a particular 
time-window, r. 

To study the market behaviour when eigenvalues take extreme values, we need 
to partition the later depending on their size. This is achieved by first normalising 
each eigenvalue in time using 

-m - (8) 

where A(r) and ct a M are the mean and standard deviation of the eigenvalues over 
a particular reference period, r. By expressing the eigenvalue time-series in terms 
of standard deviation units (SDU), we can then easily partition each according to 
it's relative size. 
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Combining the techniques of multivariate correlation analysis and multiscaling, 
we are able to provide some novel insight into the dependence of correlation on both 
time and scale. To enable this analysis, we first show how correlations are calculated 
using wavelet coefficients. 



2.3. Wavelet covariance and correlation 

The wavelet covariance between functions f(t) and g(t) is defined to be the co- 
variance of the wavelet coefficients at a given scale. The unbiased estimator of the 
wavelet covariance at the j th scale is given by 

^ t=Lj-l 

where all the wavelet coefficients affected by the boundary are removed [33] and 
Mj = N—Lj + 1, The wavelet variance at a particular scale Vffe) is found similarly. 

The MODWT estimate of the wavelet cross correlation between functions f(t) 
and g(t) may then be calculated using the wavelet covariance and the square root 
of the wavelet variance of the functions at each scale j. The MODWT estimator 
[34] of the wavelet correlation is then given by: 

= Z 7rtt \ (10) 



3. Data 

The initial data set comprises the 49 equities of the Dow Jones (DJ) Euro Stoxx 50 
where full price data is available from May 1999 to August 2007, resulting in 2183 
daily returns, Fig. [T{a). The Dow Jones Euro Stoxx 50 is a stock index of Eurozone 
equities, designed to provide a blue-chip representation of Eurozone supersector 
leaders. The small number of stocks in the index allows calculation of the cross- 
correlation matrices for small time windows, without reducing the matrix rank. 

The second data set studied consists of high-frequency returns for the Dow Jones 
(DJ) Euro Stoxx 50 from May 2008 to April 2009, resulting in 109, 540 one-minute 
Returns, Fig. [Ub). High-frequency data allows us to examine additional features 
and analyse whether the correlation dynamics previously found, [llj . are inherent 
at all time-scales or if there is a gradual emergence. The time-frame studied is 
of particular interest, with sustained market drops and high levels of volatility. 
Understanding market interactions in a period such as this may assist in forming 
portfolios that are robust to large downwards moves. 



4. Results 

To decompose the data into component time-scales, we selected the least asymmet- 
ric (LA) wavelet, (known as the Symmlet), which exhibits near symmetry about 
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DJ Euro Stoxx 50, Daily Prices, May 1999 - June 2008 



DJ Euro Stoxx 50, 1 Minute Prices, May 2008 - April 2009 




Fig. 1. a) DJ Euro Stoxx 50, Daily Prices May 1999 - June 2008 (b) DJ Euro Stoxx 50, 1 Minute 
Prices May 2008 - April 2009 



the filter midpoint. LA niters are available in even widths and the optimal filter 
width is dependent on the characteristics of the signal and the length of the data. 
The filter width chosen for this study was the LA8, (where the 8 refers to the width 
of the scaling function) since it enabled us to accurately calculate wavelet correla- 
tions to relatively long scales, dependent upon the time window used. Although the 
MODWT can accommodate any level Jo, in practise the largest level is chosen so as 
to prevent decomposition at scales longer than the total length of the data series. 



4.1. Medium Frequency Eigenvalue dynamics 

We first focus on medium-frequency equity returns measured at a one-day hori- 
zon. For each stock, using a sliding time-window of 100 days, we decompose the 
returns into their component scales in each window, using the wavelet transform. 
The correlation between the wavelet coefficients, corresponding to each equity, is 
then calculated at each scale, as described (Section HO)) . Time-series of eigenvalues 
are found at each scale, by calculating the associated eigenspectrum. This approach 
of calculating the eigenspectrum in sliding time-windows, allows analysis of the 
dependence of correlation dynamics on granularity. 

The time-window was chosen such that Q = £ = 2.04, thus ensuring that the 
data would be close to stationary in each sliding window, (different values of Q 
were examined, [TT], with little variation in the results). Fig. [TJa) shows the value 
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Fig. 2. a) Largest eigenvalue dynamics original data (b) 3 day scale (c) 6 day scale (d) 11 day scale 



of the DJ Euro Stoxx Index over the period studied. Fig. [SJa) displays the largest 
eigenvalue, calculated using the unfiltered (one day) time-series data for sliding time 
windows. As shown, the largest eigenvalue is far from static, rising from a minimum 
of 7.5 to a maximum of 30.5 from early 2001 to late 2003, (coinciding with the 
bursting of the "tech" bubble). This corresponds to an increase in the influence of 
the "Market" , with the behaviour of traders becoming more correlated. The next 
major increase occurred in early 2006, followed by a relatively marked decline until 
the beginning of the "Credit Crunch" in 2007. Similar to [17] , we note an increase in 
the value of the largest eigenvalue during times of market stress, with lower values 
during more "normal" periods. 

We next calculate, using the MODWT (Section 12. ip . the value of the largest 
eigenvalue of the cross-correlation matrix over longer time horizons of 3, 6 and 11 
days (Fig.|2Jb-d)). Certain traders, (such as Hedge Fund managers), may have very 
short trading horizons while others, (such as Pension Fund managers), have much 
longer horizons. By looking at the value of the largest eigenvalue at different scales, 
we try to characterise the impact of these different trading horizons on the cross- 
correlation dynamics between large capitalisation stocks. In Fig.[^b-d), we see that 
the main features found in the unfiltered data are preserved over longer time scales. 
However, certain features, such as the sizeable drop in the largest eigenvalue at the 
longest scale in late 2003, are not seen at shorter scales, but the aggregate impact for 
the unfiltered data is a moderate drop. Other features, such as the increase in 2006, 
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are not preserved across all scales. The different features, found at various scales, 
suggest that the correlation matrix is made up of interactions between stocks, traded 
by investors with different time horizons. This has implications for risk management, 
as the correlation matrix used for input in a portfolio optimisation should depend 
on the investor's time horizon. We look more closely at this in Section 14751 

4.2. High-Frequency Correlation Analysis 

Using the entire high-frequency data set for the Euro Stoxx 50, described earlier, we 
calculated the average correlation and eigenvalue spectrum at a number of scales, 
Table [T] As expected, the largest eigenvalue increases as the average correlation 
between stocks increases, in keeping with the 'one-factor' toy-model of correlations 
described previously, [llj . Analysis of Table [1] gives some intuition for the results 
shown, Fig. [21 with increases in the largest eigenvalue in time corresponding to an 
increase in the average correlation between stocks, ('single factor' model, |llj). 





Time 


Average 


■yst 


2 nd 


3 rd 


Scale 


Horizon 


Correlation 


Eigenvalue 


Eigenvalue 


Eigenvalue 


1 


1 min 


0.35 


20.24 


6.66 


1.085 


2 


5 min 


0.39 


20.66 


5.53 


1.26 


3 


30 min 


0.45 


23.33 


4.229 


1.64 


5 


120 min 


0.47 


24.7 


4.18 


1.74 


8 


480 min 


0.52 


27.72 


3.54 


2.42 


10 


1100 min 


0.51 


27.15 


3.83 


1.97 


12 


2000 min 


0.41 


18.42 


2.78 


2.39 



Table 1. Correlation and Eigenspectrum Analysis. The 1 st eigenvalue is the largest eigenvalue of 
the correlation matrix of Euro Stoxx 50 high frequency data from May 2008 to April 2009. 



As mentioned, (Section [1}, the Epps effect is the gradual build-up of correlation 
as the scale at which stock returns are measured, increases. In previous work, [291 
[30], the Epps effect has been studied using either the bivariate correlation between 
pairs of stocks or the average of a number of such pairs. However, as seen in Table[TJ 
this effect is also visible using a multivariate approach, with the largest eigenvalue 
increasing from 20.24 at the smallest scale to 27.72 at the 480 min (w 1 day) scale, 
with corresponding increase in average system correlation (from 0.35 to 0.52). 

In previous studies, Authors have calculated the correlation at different scales 
by recalculating the returns, resulting in reduced quantities of available data, and 
constraining the maximum possible scale. However, using Wavelet Multiscaling to 
decompose the data into component scales, we can calculate the correlations and 
eigenvalues at much longer scales. Of particular interest is the increase in the av- 
erage correlation to a maximum at a scale of 480 minutes and the subsequent 
decrease at longer scales of 2000 min (« 4 days) to 0.41. This decrease in correla- 
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tion corresponds to a decrease in the 'Market Effect', with traders having longer 
time-horizons acting in a less synchronous fashion. The decrease may be due to the 
additional discontinuities incorporated in the data at scales longer than one day as 
additional information is incorporated as stock markets open for the day, (visible 
in discontinuous jumps in the equity prices). 

4.3. High- Frequency Correlation Dynamics 

Building on the analysis of previous Sections, we now examine the dynamics of the 
eigenvalues over time using the high-frequency data described. The methodology 
is the same as for daily data, described in Section 14.11 using a time window of 
10,000 minutes. The results for small scales are shown, Fig. 03a), with a general 
increase in the largest eigenvalue for longer scales, corresponding to the increase in 
correlation shown in Section |4~2"1 However, for small scales (1 — 5 minutes), some 
distinct dynamics are found, with a sharp increase at around 50, 000 minutes and 
a decrease at longer scales (30 — 60 minutes). These distinct dynamics emerge even 
more distinctly at longer scales, Fig.(3](b), with eigenvalues actually having opposite 
behaviour at certain points, (the largest eigenvalue increases at the 1 minute scale 
but actually decreases at the longest scales, at approximately 50,000 minutes). 



Euro Stoxx 50 Correlation Dynamics, Short Scales 



(a) 




01 23456789 10 

Time (Minutes) x ^q 4 



Euro Stoxx 50 Correlation Dynamics, Long Scales 



(b) 




0123456789 10 

Time (Minutes) x10 4 



Fig. 3. Eigenvalues dynamics 

Of particular interest is the gradual increase in the largest eigenvalue moving to 
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longer scales, in keeping with the Epps effect. As found in the previous Section at 
longer scales, (2,000 minutes), using the complete data set, the largest eigenvalue 
is smaller than for shorter scales at certain points in time. This variation in the 
largest eigenvalue corresponds to a change in the average correlation between stocks. 
The distinct behaviour at longer time scales may be due to a number of factors, 
with longer scales incorporating additional information and smoothing out high 
frequency noise. 



Euro Stoxx 50 Correlation Dynamics, 2 nd Largest Eigenvalue 




0123456789 10 

Time (Minutes) x 10 * 




Fig. 4. Eigenvalues dynamics, 2 nd and 3 rd largest eigenvalues 



The dynamics of the 2 nd and 3 rd largest eigenvalue are also shown in Fig. 01 
Again, the relative size and dynamics of the eigenvalues are found to be scale- 
dependent, with a broad increase in the eigenvalues for longer scales. This is in sharp 
contrast to the results found in Section |4~21 where the 2 nd largest eigenvalue was 
found to strictly decrease as the scale increased. It appears that for small windows, 
the second largest eigenvalue is of greater relative importance (at longer scales) with 
a corresponding redistribution of correlation structure across the eigenspectrum, 
(eigenvalue repulsion). In the context of previous work done on Random Matrix 
Theory, [H[§], this may imply that different information is to be found in the 2 nd and 
subsequent eigenvalues at longer scales. These eigenvalues are said to correspond to 
different sectoral groupings, with possible group changes as we move down the scales. 
By removing the market mode, researchers [24| . found the correlation structure to be 
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invariant across scales. However, this paper did not examine different time- windows 
or longer time-scales (> 255 Mins). By combining the method of these authors with 
those proposed above, it may be possible to investigate the invariance of correlations 
across time as well as scale. 

4.4. Drawdown analysis 

As indicated, (Section [T] and [TTJ |T7]), drawdowns, (or periods of large negative 
returns), and drawups, (periods of large positive returns), tend to be accompanied 
by an increase in different eigenstates of the cross-correlation matrix. In this sec- 
tion, we attempt to characterise the market according to the relative size of the 
eigenvalues, as well as through the use of eigenvalue ratios. 

The returns, correlation matrix and eigenvalue spectrum time-series for over- 
lapping windows of 10, 000 minutes were calculated and normalised using the mean 
and standard deviation over the entire series, (Eqn. [8]). By representing normalised 
eigenvalues in terms of standard deviation units (SDU), we can partition the eigen- 
values according to their magnitude. The average return of the index is shown in 
TableEl during periods when the largest eigenvalue is ±1 SDU). As the market trend 
was predominantly negative, average returns were negative both for periods when 
the largest eigenvalue was equal to its mimimum and maximum values. However, 
distinct behaviour still emerged across all scales. 



Scale 


No. Std > 1 


< -1 


Original 


-8.06% 


-3.39% 


5 Minutes 


-7.95% 


-3.06% 


60 Minutes 


-9.27% 


-3.64% 


360 Minutes 


-8.41% 


-3.52% 


1100 Minutes 


-10.79% 


0.01% 


2000 Minutes 


-12.06% 


-6.32% 



Table 2. Drawdown/Drawup analysis. Average Index Returns when various eigenvalue partitions 
in SDU are > 1 and < — 1. 

Looking first at the original or unfiltered data, we see that when the largest 
eigenvalue is > 1 SDU, the return of the index was —8.06%. In contrast, when the 
largest eigenvalue is < —1 SDU, the return was —3.06%. Moving up the scales, 
similar behaviour was found with large eigenvalues being characterised by more 
pronounced downward moves. The emergence of this characteristic at very small 
scales (« 5 minutes), implies that it is not a result of data granularity. 

4.5. Portfolio Optimisation 

As mentioned previously, the time-scale dependence of correlation structure has 
implications for risk management. In this Section we demonstrate, using high- 
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frequency 1 minute stock returns, that the optimal portfolio found from classic port- 
folio optimisation is time-scale dependent. The correlation matrix between stocks 
of the Euro Stoxx 50 was calculated for two time-frames, the first being a short 
time window of 10, 000 minutes, (corresponding to the first time window studied in 
Section l4~3j) . the second the complete data set of 109,540 minutes, (corresponding 
to the data studied in Section |4.2|) . The portfolio optimisation was unconstrained, 
allowing 'short-selling' of stocks, which means that positive expected returns for the 
portfolio were possible, even though the majority of expected returns were negative. 
The expected returns and risk for the equities were the same for both optimisations, 
meaning that all deviations in the risk/return profile are due to changes in the cor- 
relation structure. 



Efficient Frontier Different Scales Using Time Window of 1 0,000 Minutes Efficient Frontier Different Scales Using Time Window of 1 09,000 Minutes 

4.5 | , , , , , 1 4.5 1 1 1 , , , , , 1 



(a) (b) 




10 20 30 40 50 60 10 20 30 40 50 60 70 80 

Risk (%) Risk (%) 



Fig. 5. Portfolio Optimisation 



The portfolio optimisation was performed in a standard fashion, [351 1371 138] . 
and the results, shown in Fig. [5l are highly dependent on the time-scale studied. 
Interestingly, for the smaller time-window of 10, 000 minutes, the portfolios tended 
to be less risky at longer scales. In contrast, using the larger time-window to cal- 
culate the correlations between stocks, the least risky portfolio was found using a 1 
minute scale. Portfolios calculated at the 60 — 750 minute scales were the riskiest, 
with longer-period scales found to be less risky. This is in keeping with the find- 
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ings of Section 321 where the average correlation peaked at a scale of 480 minutes, 
meaning that further diversification benefits (lower correlation) can be found at the 
shortest and longest scales. 

The optimal portfolio is clearly dependent on the time-scale used to calculate the 
correlation matrix, which in turn, is dependent on the time window used. Obviously, 
this means that care must be taken by investors in the choice of both time- window 
and granularity used in the calculation of correlation matrices. 



5. Conclusions 

Using a multivariate technique first applied to complex EEG seizure data, the corre- 
lation structure between medium and high-frequency financial time-series was stud- 
ied. Application of the Maximum Overlap Discrete Wavelet Transform extended the 
study of changes in the cross-correlation structure across both scale and time. 

(1) Using the MODWT and a sliding window, the dynamics of the largest eigen- 
value of the correlation matrix were examined and shown to be time dependent 
at all scales, (using both medium and high-frequency equity returns). Similar 
dynamics were visible across all scales, but with particular features markedly 
apparent at certain scales. This suggests that the correlation matrix between 
equities consists of interactions caused by traders with different time horizons. 

(2) Study of the 2 nd and 3 rd largest eigenvalues over both time and scale revealed 
a large increase in sectoral correlations for longer scales, using high-frequency 
data. 

(3) High frequency data were used to calculate the average correlation and largest 
eigenvalue for various scales. The Epps effect was found using this multivariate 
technique and analysis of longer scales revealed a drop in the average system 
correlation at scales beyond one day. 

(4) A partition of the time-normalised eigenvalues demonstrated quantitatively the 
relationship between the size of the largest eigenvalue and the return of the 
index. This feature emerged at very small scales (< 5 minutes), implying that 
this relationship is not due to the granularity of the data studied. 

(5) Finally, we looked at the problem of portfolio optimisation at various scales. 
Results show that the optimal portfolio depends not only on scale but on the 
time-window used in the calculation of the correlation matrix. 

The correlation technique used in this paper to measure the interaction between 
agents suffers from the drawback of being linear and hence neglects any higher-order 
relationships. Future work includes the application of non-linear information-theory 
based dependence measures, which will allow the detection of complex changes in 
synchronisation behaviour around extreme financial events. These non-linear tech- 
niques may uncover further emergent features to those highlighted above and may 
result in signals that warn of likely future market turmoil. The interaction of various 
agents with competing strategies and operating over different time scales, requires 
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further investigation in order to improve understanding of complex phenomema 
such as financial crashes. Additionally, while this financial work has benefited from 
parallels drawn with EEG system analysis, insights gained here may also raise ques- 
tions with regards to the behaviour of other complex systems. 
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